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Abstract. Recent results obtained within the Hamiltonian approach to continuum Yang-Mills the- 
ory in Coulomb gauge are reviewed. 



INTRODUCTION 

There have been many attempts in the past to solve the Yang-Mills Schrodinger equation 
using gauge invariant wave functionals. Unfortunately, all these approaches have not 
been blessed with much success 2 . The reason is that it is extremely difficult to work with 
gauge invariant wave functionals. A much more economic way is to explicitly resolve 
Gauss' law (which ensures gauge invariance of the wave functional) by fixing the gauge. 
For this purpose, the Coulomb gauge d ■ A = is particularly convenient. In my talk, I 
would like to present a variational solution of the Yang-Mills Schrodinger equation in 
Coulomb gauge 02D. Let me start by briefly summarizing the essential ingredients of the 
quantization of Yang-Mills theory in Coulomb gauge. 

In Coulomb gauge the space of (transversal) gauge orbits has a non-trivial metric, 
which is given by the Faddeev-Popov determinant 7(A) = Det(-Didi), where t) ab = 
8 ab dj +Af b , A ab = f acb Aj denotes the covariant derivative in the adjoint representation 
of the gauge group (f acb is the structure constant). In Coulomb gauge the Yang-Mills 
Hamiltonian is given by ^ 

H= 1 -Jj- l UJU+ 1 -jB 2 + ^- J J- l p(-Dd)- l (-d 2 )(-Dd)- l Jp , (1) 

where Tlf (x) = 8/ i5Af(x) denotes the momentum operator, representing the color elec- 
tric field, B is the color magnetic field and p a (x) = —Af (x)U b (x) is the non-Abelian 
color charge of the gauge field. The first term in the Hamiltonian is the Laplacian in 
curved space, the second term represents the potential and the last term arises from the 
longitudinal momentum part of the kinetic energy after resolving Gauss' law. This term 
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is usually referred to as the Coulomb term, since in the Abelian case it reduces to the 
ordinary Coulomb potential. 



VARIATIONAL SOLUTION OF THE YANG-MILLS 
SCHRODINGER EQUATION 

We have performed a variational solution of the Yang-Mills Schrodinger equation for 
the vacuum using the following ansatz for the wave functional J2] 

¥[A] = ^K/[A] _a exp \-\J A(oA ^j ( 2 ) 

with 3 a = i, where 0)(\x — is the variational kernel determined by minimizing the 
energy 

( V F|#| V I') = j @AJ[A]W*[A]HW[A] . (3) 
This gives rise to a set of coupled Dyson-Schwinger equations for the ghost propagator 

G ab {x,x') = (¥\(x,a\(-fid)- 1 \x',b)\ x ¥) = -{-d 2 )- l d{x,x') (4) 

g 

with d being the ghost form factor, and the gluon propagator 

D${x,x') = (y\A?(x)A b j(x)\y) = isOttjWa-^xJ) , tijix) = 8ij -X^-. (5) 

The gluon Dyson-Schwinger equation has the form of the dispersion relation of a 
relativistic particle 

co 2 (k) =k 2 + Z 2 (k) , (6) 
where the gluon self-energy E is dominated by the ghost loop 

which is a measure for the curvature of the space of transversal gauge orbits and which 
in the following will be referred to as curvature. Due to the ansatz for the ghost form 
factor d the coupling constant g drops out from the Dyson-Schwinger equations which 
is a specific feature of the one-loop approximation. 



This ansatz is motivated by the form of the wave function of a point particle in a spherical symmetric 

potential in a zero angular momentum (s-)state = ^y-, where J(r) — r 2 . In ref. \4] the ansatz (|2ji with 
a = was used. 



In ref. |0] an infrared analysis of the Dyson-Schwinger equations was performed without 
resorting to the angular approximation. Using the infrared ansatze 



1 r~i 
and implementing the horizon condition d~ (0) = 0, first given in ref. [6], one finds 

from the ghost Dyson-Schwinger equation the sum rule 

a + 20=d-4, (9) 

where d is the number of spatial dimensions, see also ref. ED. The sum rule © is due to 
the non-renormalization of the ghosHduon vertex which is a feature of both the Landau 
gauge Hi and the Coulomb gauge @]. Inserting this relation into the gluon Dyson- 
Schwinger equation, one obtains for the infrared exponent a unique solution = 0.4 in 
d = 2 and two solutions in d = 3, 

01=0.796,02=1.0. (10) 

Only 02 was previously found in the angular approximation |0], while 0i corresponds 
to the numerical solutions reported in ref. |0], where =0.85 was found. Recently the 
full (numerical) solution corresponding to the infrared exponent 02 was also found iflpll 
and is shown in fig. 1. 
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FIGURE 1. The ghost form factor d(k) (left panel) and the gluon energy <X)(k) and the curvature %(k) 
(right panel) 

The infrared exponent = 1 together with the sum rule © implies an infrared 
linearly diverging gluon energy, see fig. 1, which is a manifestation of confinement, 
and also produces a linear rising static color potential, shown in fig. 2. Let me stress 
that it is absolutely crucial to keep fully the Faddeev-Popov determinant J [A] in both the 
Hamiltonian and the integration measure. The Faddeev-Popov determinant converts the 
flat integration measure / £>A into the Haar measure of the gauge group 111 ill . When 
the Faddeev-Popov determinant is neglected [4] the gluon energy becomes infrared 
finite and the linear rise in the static color potential is lost. On the other hand, the 
obtained infrared behavior is quite robust with respect to changes of the variational wave 
functional. Extending the variational ansatz, eq. 0, to arbitrary a (where U = \ and 



25 




X/X„ p/ 0)1(1) 

FIGURE 2. Left panel: The static Coulomb potential. Right panel: The three-gluon vertex at the 
symmetric point, full line with non-perturbative propagators, dashed line with perturbative propagators. 



a = correspond to the wave functional used in refs. [2] and [4], respectively), it can 
be shown that up to two loops in the energy, the variational solution and in particular the 
ghost and gluon propagators are actually independent of a, see ref. JT2] . Other quantities 
like the three-gluon vertex are, however, more sensitive than the energy and do depend 
on a. 

The leading contribution to the three-gluon vertex is given by the triangle diagram 
with an internal ghost loop. In this order the three-gluon vertex has the same color 
structure as the bare one and can be expressed by a single form factor 

r^ k c (p uP2 ,P3)=f ahc {i(pi)j5 lk F(p u p2,p 3 ) + permutations) . (11) 

The form factor F is shown in fig. 2 for the symmetric point pi = P2 = P3- The infrared 
analysis performed in ref. |0] shows F{p 2 ) ~ (p 2 ) i^ 2 ^i^ . Inserting here d = 3 and 
/3 =0.85 from ref. [2], one finds F(p 2 ) ~ (p 2 ) ~ L775 ; which is in perfect agreement with 
the numerical result shown in fig. 2, from which one extracts F(p 2 ) ~ (p 2 )~ L77 . We 
have also investigated the ghost-gluon vertex and found that, like in Landau gauge, it 
becomes bare when one of the external ghost momenta vanishes, so that its dressing can 
be ignored, at least in the infrared. 

We can define the RG-invariant running coupling constant, 

a(k) = ^p 5 G 2 (p)D(p), (12) 

from the ghost-gluon vertex 0], where G(p) and D(p) are the ghost and gluon prop- 
agators, respectively. Due to the sum rule © for the infrared exponents this coupling 
has an infrared fixed point which for the DSE solutions with jB = 1.0 is given by 
a (0) = \6n/ (37Y C ) ~ 16.76//V C , see ref. B5J. This value obtained in the infrared analysis 
is in excellent agreement with the value of a(0) = 16.94 /N c , see II 1 Oil . which is obtained 
from the numerical solution of the DSE presented in fig. 1 . Let us also stress that this re- 
sult disagrees with the value for the infrared fixed point obtained in the Coulomb gauge 
limit of the interpolating gauges |J. 



THE 'T HOOFT LOOP 



The class © of wave functionals yield (up to two loops in the energy) the same Dyson- 
Schwinger equations independent of a. To test our wave functional and to get more 
detailed information on the structure of the Yang-Mills vacuum, observables more sen- 
sitive than the energy density should be calculated. Moreover, the confining Coulomb 
potential obtained above does not necessarily guarantee that our wave functional de- 
scribes indeed a confining vacuum state. The reason is that the Coulomb potential is 
only an upper-bound to the true static quark potential ifoll . Thus a confining Coulomb 
potential is a necessary but not a sufficient condition for a confining vacuum. 

In pure Yang-Mills theory the order parameter of confinement is the temporal Wil- 
son loop, which obeys an area law in the confined phase and a perimeter law in the 
deconfined phase. Unfortunately, this quantity is difficult to calculate in the continuum 
theory due to the path ordering. Fortunately, there exists another (dis-)order parameter 
of Yang-Mills theory which is easier to calculate in the continuum theory. This is the 
spatial 't Hooft loop which is dual to the temporal Wilson loop in the sense that it obeys 
a perimeter law in the confined phase and an area law in the deconfined phase ifLUl . 
The operator V (C) of 't Hooft's disorder parameter is defined for a spatial loop C by 



where W(C) is the spatial Wilson loop, Z is a (non-trivial) center element of the gauge 
group (Z = — 1 for SU(2)) and L(Ci,C2) denotes the Gaussian linking number between 
the two spatial loops C\ and C2. Unfortunately 't Hooft did not give an explicit realization 
of the operator V(C) in the continuum theory, but defined the operator implicitly as 
having the effect of a singular (on the loop C) gauge transformation. 

The physical meaning of the 't Hooft loop defined by equation (fT3l) is recognized by 
noticing that a time independent (spatial) center vortex srf {C\ ) whose magnetic flux is 
localized on the closed loop C\, produces the Wilson loop W \s& (Ci)l (C2) = Z L ( CuC2 \ 
It follows that the 't Hooft loop operator V(C) can be interpreted as a center vortex 
creation operator. The operator V (C) which generates a center vortex gauge field srf (C) 
i.e., y(C)¥(A) = ¥(A + ^(C))js given by 



the following relation ifLHl : 





(13) 



c 




(14) 




(15) 



with £ = £ a T a being a co-weight vector defined by exp(?'£) = Z actually satisfies the 
defining relation dT3l . Since 11? (x) is the operator of the electric field, the 't Hooft loop 
measures the electric flux (projected on the co-weight £ = C, a T a ) through the closed 



loop C, while the (spatial) Wilson loop W(C) measures the magnetic flux through C. 
The operator V(C) is not manifestly gauge invariant, but produces a gauge invariant 
result, when acting on gauge invariant states ^(A^) = *P(A). 

The wave functional in Coulomb gauge satisfies Gauss' law (which ensures gauge 
invariance) and has, hence, to be considered as the restriction of the gauge invariant 
wave functional to transverse gauge fields. Therefore by implementing the Coulomb 
gauge by the standard Faddeev-Popov method, we obtain for the 't Hooft loop 

(y(C)) = J &AJ(A) } ¥*(A) } ¥{A + ^(C)) , (16) 

where srf^{C = <9E) is the transverse part of the center vortex field stf (E) (TT3T) given by 
id 

£/ ± (dZ 1 x) = -i; I da ik dl , dd ik = Z M dx e (17) 

Jdz 4n\x — x\ 

which manifestly depends only on the loop C = <9E but not on the enclosed surface 
E as the field (1131 does. It is straightforward to calculate from eq. (ITvb 't Hooft's 
disorder parameter inserting for *P(A) the wave functional determined previously J2]. 
After straightforward calculations, one finds 

/-•oo 

(7(C)) = exp(-S) , 5(C) = / dqK(q)h(C,q) , (18) 

Jo 

where 

^^X'-S!) (19) 

contains the whole information about Yang-Mills vacuum while the quantity h(C,q) 
is exclusively determined by the geometry of the considered loop. To simplify the 
calculations of this quantity, we consider a planar circular loop with radius R. A detailed 
analysis shows that the R dependence of S(R) is determined by the infrared behavior of 
K(q). Given that G)(q) and %(q) are both infrared divergent and differ only by a finite 
constant, one finds 

K(q - 0) = (C0(q) -X{q)) q .+o = c ( 2 °) 

where c is a finite renormalization constant which can be chosen arbitrarily. Choosing 
c = 0, which implies that the infrared wave functional reduces to the strong coupling 
limit *P[A] = 1, one finds indeed a perimeter law S(R) = kR, while for c ^ one obtains 
Si? ~ R \\\(R/Rq). Finally, if one ignores the curvature (% = 0) in eq. (fT9l) . one finds an 
area law S(R) = pR 2 . Thus again we find the curvature to be crucial for the confinement 
properties of the theory. 



SUMMARY AND CONCLUSIONS 



In my talk I have discussed recent results obtained in the Hamilton approach to Yang- 
Mills theory in Coulomb gauge. Using ghost dominance, an infrared analysis of the 
Schwinger-Dyson equations has revealed both quark and gluon confinement. The gluon 
energy is infrared diverging implying the absence of asymptotic gluon states from the 
physical spectrum, and the static Coulomb potential is linearly rising. These results are 
confirmed by a full numerical solution of the Schwinger-Dyson equations. While the 
ghost-gluon vertex remains basically unrenormalized and in particular becomes the bare 
one when the ghost momentum vanishes, the three-gluon vertex is strongly infrared 
divergent. As a first non-trivial test of the variationally determined vacuum wave func- 
tional, we have calculated 't Hooft's disorder parameter choosing the renormalization 
constants of the gap equation in such a way that the wave functional reproduces in the 
infrared the strong coupling limit. Then the 't Hooft loop shows indeed a perimeter law, 
indicating a confining vacuum. 
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